Coherent lunar effect on solar neutrino:double slitt interference 
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Coherent interaction of solar neutrino with the moon and its implication are investigated. We 
solve the Dirac equation with the moon potential and show that a phase shift of the neutrino 
wave function becomes almost one unit if the neutrino penetrates the moon. Spatial interference 
is generated for the neutrino wave packets during eclipse and possibility of observing interference 
effect that is similar to that of a double slit experiment is pointed out. 
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Since neutrino from Supernova were observed [1], neutrino astronomy started. An astronomical observation using 
neutrino could supply a new information which is not obtained by light. However neutrino interacts with matter so 
weakly that event rate is extremely small. High luminosity and enhancement of event rate are required to make a 
real observation possible. The most luminous source of neutrino around the earth is sun and it is practical to use 
O ,■ solar neutrino for astronomical purpose. Especially solar neutrino of discrete energy from Be process can be a useful 
source for astronomical observations. Coherent scattering of neutrino by neutral current interactions can enhance the 
■ event rate. We study a coherent moon effect to solar neutrino in this paper. 
t-H ' Coherent interaction of neutrino with massive objects is a unique macroscopic quantum phenenomenon which 
appears despite an extreme weakness of interaction with matters. The flux of solar neutrino is 10 9 / cm 2 sec for Be 7 
neutrino and in a charged current event, which occurs incoherently, a rate is about one hundred per year in Super 
Kamiokande. A mean free path of 1 MeV neutrino in the earth is about 4 x 10 10 Km and is slightly smaller in 
| the moon. The earth and the moon are quite transparent as far as the charged current interaction is concerned. 
. Strength of the neutral current interaction is smaller than the charged current interaction and elastic scattering 
process v e + e~ — > e~ + v P is similar to the neutral current interaction but has a larger strength. In these processes 
the neutrino and other matters can stay in almost the same quantum states, and by the coherent sum of amplitudes a 
total amplitude is enhanced. This effect becomes important when neutrino interacts with a massive object like a star. 
The moon is a star around the earth and can be located between the earth and the sun. So the coherent interaction 
of the solar neutrino with the moon can take place. We study the coherent scattering of the solar neutrino with moon 
in this paper. 

The MSW coherent matter effect Rcf. [2-3] is actually one of the important ingredients in understanding current 
neutrino oscillation experiments. Because expected day-night effect has not been observed, the MSW effect would not 
be important for the earth and for the moon. Double slitt like effect we study in this paper is a spatial interference 
effect of the lightest neutrino. This effect is sensitive to the absolute magnitude of the lightest neutrino mass. 

The solar neutrino has an energy in MeV region and the coresponding de Broglie wave length is 10~ 13 m ,whearas 
the radius of the moon is of order 10 6 m. Because the wave length is much smaller than the total size, this regime 
may be thought as a geometrical optics regime or the classical regime for the neutrino where the wave character is 
normally washed out. However the interaction strength of neutrino with matter is extremely weak and the neutrino 
propagates almost freely even in matter with keeping an original phase. Consequently the propagation of neutrino in 
the moon should be treated by quantum mechanical wave equation even in the macroscopic regime. Interference is a 
key phenomenon of wave dynamics. We show that to generate interference for solar neutrino, wave packet character 
, which have been studied for neutrino oscillations in Ref. [4-8], is important. Using reasonable parameter values, 
we estimate the magnitude of the effect from overlapp between an initial wave packet and a final wave packet. We 
show that an observation of the effect would be possible if the absolute value of the neutrino mass is 10~ 4 eV or less. 
Observation of the interference should verify quantum mechanics in Gm region and give new informations on the 
neutrino mass and on the interior of the moon. 
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By replacing the quark pair operators or the electron pair operators in the four Fermi interaction with the density of 
electrons and the quarks, we obtain a potential term in Dirac equation. We solve the equation and find that the effect 
of the potential arises in the momentum dependent phase factor. Although the effect disappers from the magnitude of 
a plane wave, the potential gives an interference effect to the wave packet since the wave packet is a linear combination 
of plane waves. Interference among the plane waves is generated and an observable effect is obtained. It is predicted 
that the solar neutrino flux changes with time when they are blocked by the moon during eclipses. 

potential 

Neutrino interacts with quarks q(x) or with the electron e(x) in matter by the charged current interaction 

^(x) 7m (1 - l5 )e(x){d(x)y»(l - l5 )u(x) + e(x)>f(l - js)v{x)} t (1) 

and by the neutral current interaction, 

Gf 

^=P(x)(l - lb ) ltl v{x){gv{u)u{x)^u{x) + g v (d)d{x)^d(x) (2) 
+gv(e)e(x)~f ,1 e(x) + axialvectors}. 



Throughout this paper we assume that the electron neutrino is lightest and the mixing angle is negligible. The 
second term in charged current interction has a same final state and initial state as electron neutrino neutral current 
interaction and by Fiertz transformation the interaction is written as 



G 



v(x){l — 'ys)j lll v(x){e(x)j fJ 'e(x) + axialvector}. (3) 



The coefficents gv(u), gv(d), gy(e) are given by the Weinberg angle 6\y as 
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9v{u) = - - 2 x - x sin 2 9 w (4) 
gv(d) = -^-2 x(-^)x sin 2 9 W (5) 



9v(e) = +2 x sm 2 6 w (6) 

Using the current value sin 2 9w = 0.23108 [9], we have 

gv{u) = 0.19, g v (d) = -0.35, g v (e) = -0.034 (7) 

The charged current interaction,Eq.Q, transforms the neutrino to the electron. Since the final state in the matter 
is different from the initial state, the process occurs incoherently. The mean free path of the neutrino in the earth 
from the incoherent scattering v + n — > e + pis computed from the crossection, a, and the neutron density p n . By 
using the crosssection and the density, 



Gp 2 



2M n E u (8) 



r 



Pn =0.556-N Avo /cm 3 , (9) 
where r is the specific gravity and Na V o is Avogadoro Number, we have the mean free path, 

I = -1 

ap 

= lO^Km (10) 

for r = 5 and E v = 1 Mev. Hence the mean free path is much larger than the size of the earth, 10 5 K m or the size 
of the moon, 2 x 10 3 Km. Neutrino propages in the earth or in the moon almost freely as far as the charged current 
interactions are concerned. 
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The strengths of neutral current interactions, Eq.©, are about the same as the charged current interactions. So a 
probability of incoherent scattering event from the neutral current interaction is also about the same as the charged 
current interaction and is very small. 

Coherent scattering is possible in the neutral current interactions since the final state in the matter is the same 
as the initial state. In the coherent scattering, amplitude from each atom is coherently added and a total amplitude 
can become much larger than the original amplitude. In this case, an enhancement is possible. We will see that the 
charge of vector current can have a coherent contribution. 

To study weak matrix elements, expectation value of the current in the matter 

(Matter 1 /(0)tJ7(0) | Matter) (11) 

= E^r/*(o) 

i 

is needed where summation over atoms is taken. Expectation values of the currents in one particle state with a 
momentum p and a spin s are 

u(p,s)j u(p,s) = —, (12) 
m 

u(p,s)jiu(p,s) = — , 
m 

w(p> s)7o75u(p, s) = 0, (13) 
«(p> s)7i75«(p, s) = Si. 

Only the zeroth component of the vector charge is positive definite and nearly 1 in the forward scattering region of the 
neutrinos. Other components are not positive definite and are small. Hence after the summation over whole atoms, 
only the vector charge becomes non-zero value J2i ^(P; s )lou(p, s) = p and others vanish. 

Consequently coherent neutrino interaction with matter is obtained by the vector charge and is reduced to an 
effective two body Hamiltonian, 

v{x)^^^v(x)V (x), (14) 

Vq(x) = V2G F {p e (x) - 0.5p neutr on{x)), (15) 

where p e (x) is the density of the electrons and p n eutron( x ) is the density of the neutron . 
plane wave 

We solve a Dirac equation with the potential term that is produced by the interaction with matter. Since the 
potential is caused by the weak interaction, the magnitude is very small and is proportional to the Fermi coupling 
constant. The range of the potential we study,on the other hand, is very large. So it is interesting to see if the effect 
of the potential is observable. 

Let us solve a Dirac equation with a left-handed potential term, 

ih—ip(x) — (id ■ p + mf3)ip(x) + V(x)( — - — )ip(x), (16) 

V(x)=V ,r<R (17) 
V(x) = 0,r>R, 

here R is a large macroscopic value and Vo is a small value. We will see that the product VoR is order 1. 
We obtain a stationary solution of the energy E 

%l)(x) = exp (Et/ih)ip{2) (18) 

with a boundary condition at z — > — oo 

<t>(x) = e lis u(k), (19) 
where u(k) is a free Dirac spinor of a momentum k. ip(x) satisfies an intergral equation, 

■4>(x) = (j)(x) + f dx'D(x - x' )V(x')(^-^)tP(x'), (20) 
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with the Green function of the Dirac operator, 

v ; J (2tt) 3 E 2 -p 2 -to 2 +ie y ' 

By applying Born approximation, we have an solution 

ip(g) = 0(f) + ip(x) {1) +ip(x) {2) + ip(x) (3) + -, (22) 
V^(f) ( ' +1) = Jdx'D{x-x" )y(f')(i_21)^(f' )«, 

*(1) _ Jk-X f d P i(p-k)-x 1 



where 



^-^sr^ro (23) 

(i— ^(2E(k) + a-(p- k)) + l5 f3m)u(k)V(k -p), 



V(k-p) = / dx'e^-P>- s 'v{x') (24) 



ATTV R 3 V() (q), 
1 sin q 

q = R\(p-k)\ 



1 sing cosg 

- r )> (25) 



Identity, 



(£(fc) + a • p + m/3)- — —u(k) 



2 

was used in Ea. (|24Jl . Finally we have 



' " 75 2E{k)u{k) + j 5 mf3u(k) + ^-—^-a ■ (p - fc)u(fc) (26) 



^(D = ^ { l _Jl 2 E(k)^u{k)F(x 1 k) (27) 

Z Z7T 

+75/3m-^2-u(fc)F(f, fc) + 1 2 75 ^ a ' ft (^kMfc)}, 

F{x,k)=[dq =— v {q)e^ (28) 

J —2k ■ q + ie 

Fi(x,k) = [ dq qi ^ v (q)e i ^ . (29) 

J — 2k ■ q + ie 



Eqs. and (J2^ become 



F(£,fc) =4^-^/1-^., (30) 
Fi(x,k) = ^-F(x,k) (31) 

in the region A; • f > i?, 1 > where £* = x/R,^t = £ — k(k ■ £)/fc 2 . In the region where R and ^ are large, the first 
term is dominant over the second term and the third term in Ea. (|27|l . and we have 

^ (f )(D = ] -Jl 2 E{k)^F{x,k)4>{x). (32) 

Z Z7T 

Thus the first order term changes with position in the range R even though the wave length is microscopic size. It 
should be noted however that this term is pure imaginary. It is suggested that correction terms modify only phase 
factor of the wave function. We show in the following that this is the case in fact. 
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The second order term is computed in a similar manner. The dominant term is given as, 



TO" = 1 ' 

where the coefficient is given as 



m {2) = {^2E{k) V ^fF^{x,k)m, (33) 



F^(x,k) = I dq lC lq 2 * 1 Vo { qi )^l vviq^e 1 ^ 2 ^. (34) 

— 2k ■ (<fi + q 2 ) + it —2k ■ q 2 +ie 



By writing the integral with a symmetric manner, we have 



F (2 \x,k) (35) 



1 f ,~ J- 1 1 1 X 

— / dq\dqi ^ ( - 1 - ) 

— 2k ■ (<fi + qz) + ie — 2k ■ q 2 + ie —2k-q\+ie 



2! 



xv (qi)vo{q2)e R 

= ^7 / d Qi d Q2 Vo{qi)Vo(q2)e R 

2! J -2k-qi+ie-2k-q 2 + ie 

1 - 2 

= 2f^,*) ■ 



Thus the second order term becomes, 



H*) {2) = |(^2^(fc)^F(f, k)fm (36) 



and higher order terms are written in the same manner 

m (l) - \{ ] ^2E{k)^F{xMm- (37) 
Adding all higher order terms we have the wave function 

ip{x) = exp{ik-x + i^^2V ^Jk 2 (R 2 -x 2 ) + (k-x) 2 )u(k). (38) 

Substituing the Fermi coupling constant, the radius and the density of the moon assuming the specific gravity r, 

Gp = 1.16 x l(T 5 (GeV0- 2 (39) 
R = 1.74 x 10 3 Km (40) 

Pe = Pneutron = 0.556-N Avo / ' CTT? (41) 

to Ea. H15(l . the numerical constant in Ea . I|38|l becomes for r = 5 

2V Q R = 0.95. (42) 

This is an interesting value to see an interference effect. However the correction is in the phase factor and disappears 
in the ^(s)! 2 . So the neutrino flux is the same as free wave if the neutrino is a plane wave. It is impossible to observe 
the interference effect using any plane wave in the present situation, 
wave packet 

Realistic solar neutrino is not a plane wave but is a linear combination of plane waves, a finite wave packet Ref.[4- 
9,11-12]. The scattering amplitude is the overlapp between the initial wave packet and the final wave packet. The 
former is determined from the beam and the latter is determined from the detector. 

From the Eq. l|38|l . the phase of the initial wave function through the moon depends on several parameters such as 
momentum and position and the wave packet is certain linear combinations of these waves. We will show that the 
overlapp between these wave packets and the final wave packets have interference effect if the mass of the lightest 
neutrino is about 10 -4 eV/c 2 or less. 

A wave packet expands during a propagation [12,13]. Speed of expansion is determined by the velocity variance and 
is dominant in the transverse direction for a relativistic particle. Using the initial size Sx, the maxmimum velocity in 
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the transverse direction is given as vt — — j^g ■ The size in the transverse direction is a product of the velocity 
and the propagation time 

Ax T = v T St, (43) 

The size for the neutrino of E = O.GMeV, Sx = 10~ 10 m becomes 3 x l0 5 Km for St = 500s and 10 3 Km for St = Is. 
The former one is the packet size when it propagates from the sun to the earth and the latter is the packet size when 
it propagates from the earth to the moon. The radius of moon is 1.7 x 10 3 Km and is about the same as the packet 
size of the neutrino at the moon which is observed at the earth with a microscopic size. So the wave packet effects 
are relavant for the solar neutrino that is detected at the earth. 

The Gaussian wave packet of the variance ah at t = 0,x = Xq expands at a much later time i >> to or at a 
previous time t << —to and behaves at a distant position x as 

ij;{x,t) = Ne- l ^^ {ct)2 - {3 -^ a)2 -^ iPx -P° )2 , (44) 

P x = mc 1 =(x-X a ), 

yW-0?-^) 2 

where N is a normalization factor and po is the central value of the momentum. From the Gaussian term the wave 
packet size is determined. The phase factor is written by the use of the momentum as, 

4> = { ^Xt. (45) 

h\Pxc\ 

The phase factor F{k 1 x) is added for the neutrino which penetrates inside the moon. 

We study the overlapp of two wave packets at a certain time at around the moon. Initial wave packet which is 
emitted at a time T\ [14] in the sun and propagates toward the earth through the Moon and final wave packet which 
is detected at a time T 2 at the detector are located around the Moon at a time t = T 2 — St where St is around one 
second. The function F(k,x) is in one of the wave function and becomes smooth function of x in this region and 
computation is straightforward. 

The initial wave packet which is emitted from the sun at {Xi, T\) behaves at (x, i) as, 



i> in (cc,t) = Ne-^W-^'-V-*^-™^-^ 3 , (46) 
1 

P Xl = mc (x - Ai). 

^(t-Ti))* = 

The final wave packet detected at (X21T2) behaves at the same (x,t) as, 

Vw(x,i) = Ne- l ^^T 2 -t)y-(3~x 2 )^^(P X2 -^f+F(S)^ (47) 
Px 2 = mc (a? - X 2 )., 

^{c{T 2 -tw-{x~x 2 y 

We have chosen the situation where the wave front of the initial wave packet does not reach the Moon but the wave 
front of the final wave packet reaches the Moon. The phase in Eas. l|46|l and Q47I ) depend on the absolute value of the 
neutrino mass and momentum. For a neutrino of the mass 10 _4 ey/c 2 and the momentum lMeV/c, the phase factor 
cj) is written as 

= 5 M x IO- 3 , (48) 
x 

x = 1000-fsTm, 

and becomes negligbly small in the scale of Moon. So the phase factor of ?/>m(x) is constant 4>o and the phase factor 
of ipout(x) is negligibly small in the same scale. Gaussian factors are almost one in the same region. Then overlapp 
of the two wave functions is given as, 



(49) 



7 



Since the Gaussian terms are positive definite the integral is reduced by the effect of phase factor F(x). The reduction 
rate depends on parameters such as wave packet sizes, the energy and others. Using the Moon radius for the 
size of detecting wave packet around Moon we have the neutrino flux at about 0.85 of the normal flux value for 
Pq™ = f% ut = Px at a total solar eclipse. Reduction is almost the same for other momenta. 

We studied coherent scattering of solar neutrino by the moon and found that the moon modifies the phase of solar 
neutrino wave function. The result is surprizing from two reasons. First one is extreme weakness of potential and 
second one is the fact that the wave character is seen in extremely large scale. They became possible because, despite 
extreme weakness of potential, the volume of finite potential region is very large and the product of two values is of 
order one. The other feature is seen from the fact that the ratio between the distance and de Broglie wave length 
is around 10 20 . Normally in this regime the wave character are washed away due to a rapid oscillation of phase 
and geometrical optics regime is realized. Interference should have not been produced in this case. However in the 
present case, interference is in fact produced. This is because the time dependent phase and space dependent phase 
of the relativistic waves cancell and consequently the phase difference in large distance is not washed away and wave 
phenomena occurs. 

As an implication, we suggests that the time dependent moduration of solar neutrino flux occur during eclipse if 
the suitable detector is used. If a macroscopic wave phenomenon of neutrino is verified, this can be used for testing 
quantum mechanics in the scale of Gmcter range. New informations on the absolute value of the lightest neutrino 
mass and the interior of the moon could be supplied also. 
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